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mass squared can be extended to third order in the fermion mass. Within a 
numerical treatment we include higher Fock states up to six particles. Our 
results are consistent with all previous numerical investigations, in particu- 
lar lattice calculations. For the massive Schwinger model, we find a small 
discrepancy (<2%) in comparison with known results. Possible resolutions of 
this discrepancy are discussed. 
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I. INTRODUCTION 



The number of papers on the Schwinger |IJ] and 't Hooft [f2| M models is legion - for 
a twofold reason. On the one hand, the two models are particularly simple and become, 
for special choices of parameters, even exactly solvable. On the other hand, despite their 
simplicity, both models contain interesting physics analogous or similar to properties of 
gauge theories in higher dimensions. The models therefore serve as interesting theoretical 
laboratories for studying phenomena like confinement and chiral symmetry breaking, to 
mention only the most prominent ones. For recent reviews on the two models we refer the 
reader to [Q for the Schwinger model and to |J for the 't Hooft model and generalizations 
thereof. 

While confinement (or charge screening) is realized in the same way in both models, via 
a linearly rising Coulomb potential, the second feature, chiral symmetry breaking, is not. 
In the 't Hooft model, 1+1 dimensional QCD in the limit of large Nq, chiral symmetry is 
'almost' spontaneously broken || [[/]] and a massless bound state arises in the chiral limit 
of vanishing quark mass || ||. In the massive Schwinger model, QED in d = 1+1 || 0, 
chiral symmetry is anomalously broken, and the contribution of the anomaly to the mass 
gap survives the chiral limit yielding the free, massive boson of the massless model. This 
boson becomes interacting in the massive model and again can be viewed as a bound state 
of fermionic constituents. For the sake of brevity we will call the lowest bound state of 
both models the 'pion' (although, in the real world of d = 3+1, the Schwinger model boson 
corresponds to the //-meson). In this paper we will try to determine its mass and (light- 
cone) wave function as accurately as possible. We mention in passing that the Schwinger 
model chiral anomaly is closely related to the appearance of a vacuum #-angle parameter 
|T0| [JTTJ] which also affects the particle spectrum ||. Throughout this paper, however, we 
implicitly assume that 6 is set to zero. 

The calculation of bound state masses and wave functions for the two models at hand 
has a long history, beginning with the exact solution of the massless Schwinger model [|TJ. 
For non-vanishing 'electron' mass, m, the model is no longer exactly solvable and one has 
to resort to approximations. One of them is to assume the mass m being small and expand 
around the massless solution [g] 0. For the 'pion' mass squared one expects an expansion 
of the form 

M 2 (m) = 1 + M x m + M 2 m 2 + M 3 m 3 + 0(m 4 ) , (1.1) 

where all masses are measured in units of the basic scale fio = e/^/w, the mass of the boson 
in the massless Schwinger model, which is thus represented by the '1' in ( |1.1| ). The first 
order coefficient Mi was obtained analytically in [jllj via bosonization, 



M x = 2e 7 = 3.56215 . (1.2) 
with 7 = 0.577216 being Euler's constant. Shortly afterwards Bergknoff, using light-cone 



Hamiltonian techniques (see below), found a value |13 



Mi = 2n/V3 = 3.62760 , (1.3) 
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which differs from ( L2|) by 1.8 %. One topic of the present work will be the analysis of this 
discrepancy and an attempt to do better by refining 't Hooft's and Bergknoff 's methods. This 
is particularly important as the coefficient Mi is directly related to the vacuum structure of 
the Schwinger model by 



Mi = -4tt(O|^|O)cos0 



;i.4) 



where (0|^?/>|0) is the condensate of the massless Schwinger model JT^] [T6|] (in units of fio), 

1 „, 



(0|^|0) 



2vr 



-0.28347 



;i.5) 



and 9 the vacuum #-angle. The result ( |1.2|) thus corresponds to 9 = 0. In [[17]], the relation 
( |1.4| ) (for 9 = 0), which is a 1+1 dimensional analogue of the Gell-Mann-Oakes-Renner 
formula fi~8l , has been used to determine the condensate from the 'pion' mass squared via 



(oHMo) 



l__d_ 

47r dm 



MHm, 



m=0 



;i.e) 



Any inaccuracy in the determination of M\ thus immediately affects the value of the con- 
densate At this point it should be mentioned that within light-cone quantization there 
have been many attempts to calculate the condensate of the massless model alternatively 
by solving for its vacuum structure ||2( 



Recently the calculations of M 2 (m) have been extended to second order. Using func- 



tional integral techniques and mass perturbation theory, Adam [^TJ found the analytical 
expression 



M 2 = 47r 2 (0|^|0) 2 (Acos2fl + J B) 



;i.7) 



where A = —0.6599 and B = 1.7277 are numerical constants, given in terms of integral 
expressions containing the 'pion' propagator for m = 0. Inserting the values for A, B and 
the condensate, and setting 9 = 0, ( |1.7[ ) becomes 



M 2 = 3.3874 . 



;i.8) 



This result has been confirmed independently by Fields et al. |22j, who derived the same 
integral expressions by summing up all relevant Feynman diagrams in the bosonized theory 
using near-light-cone coordinates. 

To first order in m, the mass-squared of the 't Hooft model 'pion' has already been calcu- 
lated by 't Hooft H by solving the associated light-front bound-state equation. He derived 
this equation by projecting the covariant Bethe-Salpeter equation onto three-dimensional 
hypersurfaces of equal light-cone time, x + . Soon afterwards, the equation was rederived 
using light-cone Hamiltonian techniques ||23[| . The result for the 'pion' mass-squared is 



71 



[m) = 2—j=m + 0(m 2 ) 
v3 



;i.9) 



As the light-front bound state equations of the two models at hand differ only by an additive 



contribution due to the anomaly (see [O] El and below), it is not too surprising that 
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't Hooft's and Bergknoff 's results coincide. For the 't Hooft model, all masses are expressed 
in units of the basic scale /xjj = g 2 Ncf2TT []. It is obvious from (|1.9f) that M 2 vanishes in 
the chiral limit, m — - > 0. As explained in [Of pi], this is not in contradiction with Coleman's 



theorem [fjl| as the 'pion' is not a Goldstone boson and decouples from the 5-matrix. 

The condensate of the 't Hooft model has first been calculated by Zhitnitsky using sum 
rule techniques @, 



(0\ipip\0)/Nc = Mi = -0.28868 

4:71 



1.10) 



As expected, the condensate is proportional to Nc and can again be written in terms of M\. 
The result has been confirmed analytically |26] |57] and numerically E3|. Via ( l.lOj) , the 
numerical value 



28 for the condensate leads to a numerical estimate for M\, 



M x = 3.64 ± 0.05 



Though the numerical calculation did not use light-front methods, it seems to favour the 
't Hooft-Bergknoff value ( |1.3| ), which is the standard value for the 't Hooft model. Higher 
order corrections to ( |I~9| ) have been discussed in |29|] without explicit calculation of the 
expansion coefficients. 

In recent years, both models have been serving as a testing ground for new techniques 
developed in order to solve bound state problems using light-cone (or, equivalently, light- 
front) quantization. These new methods are 'discretized light-cone quantization' (DLCQ) 
IPI [pi]] , where one works in a finite volume leading to an equally-spaced momentum grid, 
and the 'light-front Tamm-Dancoff (LFTD) approximation' |[32|| , which (drastically) trun- 
cates the Fock space of the continuum theory thus limiting the number of constituents a 
bound state can have. The latter approach can be viewed as a generalization of the tech- 
niques of 't Hooft and Bergknoff. Both methods aim at a numerical solution of relativistic 



bound state problems. DLCQ has been applied to the massive Schwinger model [|33|] and to 
QCD in 1+1 dimensions (for arbitrary N c ) |34] 
for both models, QCDi + i |35j |36| and the Schwinger model [j3?H 



There are analogous LFTD calculations 
These latter works 



are closer in spirit to ours than the DLCQ approaches. We will compare to all this recent 
work in more detail later on. 

The purpose of this paper is to obtain the 'pion' mass squared of both the Schwinger 
and 't Hooft model to high accuracy including the third order in m. In addition, we want 
to calculate the light-cone wave function of the 'pion' with high precision. We will use 
the light-front techniques of 't Hooft and Bergknoff and extensions thereof. Our starting 
point is a LFTD truncation to the two-particle sector. The low order calculations will be 
done analytically. To go beyond, computer algebraic and numerical methods will be applied 
and their convergence tested. The final goal is to shed some light on the advantages and 
limitations of this particular light front approach to bound-state equations. 



1 Some authors (including 't Hooft) use a different convention for the coupling, which amounts to 
the replacement g 2 — > 2g 2 . 
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The paper is organized as follows. In Section || we review 't Hooft's ansatz for the wave 
function yielding the lowest order solution of the bound-state equation. We compare the 
exact endpoint analysis with a variational procedure which is introduced at this point. In 
Section [HI] we extend 't Hooft's ansatz by adding more variational parameters in such a way 
that an analytic solution can still be obtained. To this end we employ computer algebraic 
methods which allow to treat up to five variational parameters. This is sufficient to achieve 
rather good convergence. In Section [TV] , these results are compared with the outcome of 
purely numerical calculations. We conclude the paper in Section |V| with some discussion of 
the presented as well as related work. 



II. 'T HOOFT'S ANSATZ 

Our starting point is the bound state equation of the 't Hooft and Schwinger model in 



the two-particle sector, which, in a unified way, can be written as (I7| |24 



M^x) = (m 2 - 1) M .vj 1 dy ^ 



x(l — x) Jo (x — y) 
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+ a / dx<p(x) . (2.1) 
Jo 

We will refer to this equation as the 't Hooft-Bergknoff equation in what follows. 4>{x) 
denotes the valence part of the 'pion' wave function, x and y the momentum fraction of one 
of the two fermions in the meson. The symbol V indicates that the integral is defined as a 
principal value [|3£fl . 

The parameter a measures the strength of the anomaly. In the 't Hooft model, a = 
(no anomaly), and in the Schwinger model a = 1 (representing the usual chiral anomaly). 
The scale parameters, as already mentioned, are given by /i 2 , = g 2 N c /2iT and /i 2 , = e 2 /n, 
respectively. M denotes the mass of the lowest lying bound state (the 'pion'). Our objective 
is to obtain (approximate, but accurate) solutions for M and 0. 

Upon multiplying ( j2.1| ) with <p{x) and integrating over x we obtain for the eigenvalue 

M 2 (m) = (m 2 - 1)£ - ^ + a 3-, (2.2) 

where we have defined the integrals 

I = ^ dx(p 2 (x) , (2.3) 

j^/Vg^ (2.4) 
JO x(l — X) 



o (x- y f 

dx(f)(x) . (2.6) 

Iq is the norm of the wave function, Ii and I2 are matrix elements of the mass and interaction 
term in the light-cone Hamiltonian [|TB|] in the state \<p). The integral J3 is the wave 
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function at the origin. An independent formula for the 'pion' mass-squared is obtained by 
integrating (|2.1|) over x. This results in the simple expression 

M 2 (m) = m 2 ^- + a , (2.7) 

^3 



with the additional integral, 



1 dx f}^ x . (2.8) 



Of course, for the exact wave functions, the right-hand sides of (|2.2j ) and (|2.?] ) have to 
coincide. As the wave functions cannot be obtained exactly, we will later use the agree- 
ment between both values for the mass-squared as a measure for the accuracy of our wave 
functions. To determine the latter we will use a particular class of variational ansatze. 

In his original work on the subject, || |, 't Hooft used the following ansatz for the 
wave function 

<j>{x) = x p {\ - xf . (2.9) 

This ansatz is symmetric in x <-> 1 — x (charge conjugation odd), and (3 is supposed to 
lie between zero and one so that the endpoint behaviour is non- analytic. As a non-trivial 
boundary condition one has the exact solution of the massless case, 

M 2 = a, and <f>(x) = 1 , i.e. (3 = . (2.10) 

The main effect of having a non- vanishing fermion mass is the vanishing of the wave functions 
at the endpoints implying a non-zero (3. This suggests the following series expansion for (3, 

(3{m) = (3im + (3 2 m 2 + (3 3 m 3 + 0(m 4 ) , (2.11) 

and for the 'pion' mass squared, 

M 2 = a + Mim + M 2 m 2 + M 3 m 3 + 0(m 4 ) , (2.12) 

in accordance with ( | . 1|) . 



A. Exact Endpoint Behaviour 



The exponent (3 in ( |2.9| ) can actually be determined exactly by studying the small-x 
behaviour of the bound state equation fl2.1|) . To this end we evaluate the principal value 
integral for x — > 0, 



V / dy 



y (i-vY 

(x-y) 2 



X: 



13-1 



V 



dz- 



-X 



0-1 



l-zf 
IT (3 cot IT (3 + 0(x) 



0(x) 



(2.13) 
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Plugging this into 
dental equation || 



H) and demanding the coefficient of x 13 1 to vanish yields the transcen 



m 2 - 1 + rrp cot %p = 



(2.14) 



which is independent of the anomaly a. Using this expression we can determine /3 either 
numerically for arbitrary m or analytically for small m by expanding 



P 



73 



m I 1 



7T 



10 



-m 



+ 0(m 4 



(2.15) 



Note that the second order coefficient, P%, is vanishing. Furthermore, for the exact (3 one 
has P1/P3 = —1/10, which we will use as a check for our numerical results later on. 

Our task is now to determine the coefficients Mj, % = 1, 2, 3 in ( |2.12| ). The ansatz 
leads to the following results for the for the 'pion' mass squared (p. 2D, 



M 2 (m,p) 




l + a 







(2/3 + 1)2 



(2.16) 



In the above, B(z\, Z2) denotes the Beta function. The relevant formulae for double integrals 
like I2 Q2.5p can be found in ]30] and [|T|, Appendix C. For a = 0, Q2.16 ) has also been 
obtained in [RH]. Let us emphasize that this result, which expresses the 'pion' mass squared 



as a function of the (exactly known) endpoint exponent /3, is only approximate as 't Hooft's 
ansatz Q2.9| ) for the wave function does not represent an exact solution of the bound-state 
equation. It is only the endpoint behaviour that is known exactly; in the intermediate-x 
region 't Hooft's ansatz is only an approximation that presumably becomes rather bad for 
large masses (non-relativistic limit) where the wave function is strongly peaked at x — 1/2. 
The accuracy of the result will be discussed extensively later on. 

In order to find M 2 to order m 3 we need to expand (|2.16|) to order /3 3 , as /3 itself is of 
order m. The result is 



Inserting /? from ( 2.15|) one finds 



AT (ml 



7T 



a + 2—=m 
\/3 



+ 



4 2 



9 



7T- 



C(3) 



a 1 



12 

7T 2 



m 



M 2 (m, p) = — + Am 2 + a + —(3 

ID O 

+ 4[tt 2 /3 - 3C(3) + a(vr 2 /12 - l)]p 2 

+ [jjvr 4 - 48C(3) + 4a(4 - 3C(3))]/3 3 
+ 0(P 4 ) . 



(2.17) 
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+ 



3v/3 



7T" 



tt 4 -48C(3) +4a(4-3C(3) 



m 



+ 0(m 4 



(2.18) 



Let us give the explicit results for the 't Hooft and Schwinger model, i.e. for a = and a 
=1, respectively. For a = 0, we have 



7T 



Mi = 2^ = 3.627599 
v/3 

M 2 = 4(2-C(3)-^ 



7T 



3.615422 



M s 



3v/I r-7r 4 - 48C(3)) = 0.125139 , 



and for a = 1, 



7T 



7T 



3 V 5 



Mi = 2^= = 3.627599 



M 2 = 4(2-C(3)-^ 



+ 1 " 



12 



3.399568 



7T 



3^ 



7T° 



0.389137 . 



(|tt 4 -48C(3))+4(4-3C(3)) 



(2.19) 
(2.20) 
(2.21) 

(2.22) 
(2.23) 

(2.24) 



We note that the anomaly does not contribute to the first order term which therefore is the 
same in both the Schwinger and 't Hooft model. This confirms the observation made in the 
introduction, the coincidence of (|L~3l ) and ( |1.9|) . As Mi is proportional to the condensate, 
cf. ( |1.6| ) and ( |1.10| ), the latter is also independent of the anomaly [ |T7| . 

For the Schwinger model (a = 1), the second order result can be compared with the 



analytical calculation of Adam, ( |1.8| ). Astonishingly, the relative difference is smaller than 
for the first order, namely 0.36 %. This is accidental, as we shall see below. 

As a cross check, we determine the 'pion' mass squared using the alternative formula 
Plugging in the ansatz ( |2.9| ), we find the fairly simple expression 



M 2 (m,(3) = Am 2 + a + 
Inserting the expansion ( |2.15|) this becomes 



2m 2 



M 



a + Mim + M 2 m 2 + M 3 m 3 + 0(m 4 



with the coefficients Mi explicitly given by 



Mi = 2n/V3 = 3.627599 , 
M 2 = 4, 

M 3 = Mx/10 = 0.03627599 



(2.25) 



(2.26) 



(2.27) 
(2.28) 
(2.29) 
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In ( |2.26| ) the anomaly contributes only in zeroth order (as (3 is independent of a), so that 
the Mi are the same for a = and a = 1. Comparing with ( p,19|) and (p. 22 ), we see that 
the first order coefficients coincide for the two alternative mass formulae. In addition, one 
finds that M 2 roughly coincides with the M 2 of both the 't Hooft and Schwinger model to 
within 10-20 %. The value for M3 is smaller by approximately an order of magnitude. 

We conclude that 't Hooft's ansatz works well to lowest non-trivial order in m but needs 
improvement if one wants to go further. As a preparative step for such a development 
we now introduce a variational method that can easily be extended to accurately calculate 
higher orders in m. 



B. Variational Approach 



At variance with the above we are now going to regard (3 as a variational parameter to 
be determined by minimizing the function M 2 (m,/3) of ( |2.17| ). To this end we insert the 
expansion ( [2.1 1|) into ( P-lTf ) and obtain 



M 2 (m) = a + 



+ 



+ 4 



+ 



¥ -3«3) + a(- 



1 W 



m 



8/5i/3 2 (y-3C(3))+/? 1 3 (^ 4 -48 

+ a(2(5 x (3^ - 4) + 4/^(4 - 3C(3) 
+ 0(m 4 ) . 



m 



^2.30) 



Note that to order m only the leading coefficient j3\ contributes. We will furthermore shortly 
see that the dependence of M 2 on /3 2 and of M 3 on j3 3 is only apparent. 

Solving the minimization equation, dM 2 /d(3 = 0, for the coefficients leads to 



p x = v^/vr = 0.55133 , 

(3 2 = 0.11690 + 0.065612a , 

p 3 = 0.0049077 - 0.050811a + 0.019521a 2 



(2.31) 
(2.32) 
(2.33) 



Comparing with ( [2.15|) we note that the coefficient (3\ is exact! Plugging it into (|2.30|) 
we verify the statement above that M 2 is independent of (3 2 and M 3 independent of (3%. 
Therefore, the expressions ( |2.30| ) and ( |2.18| ) coincide up to and including order m 2 . Mi 
and M 2 are thus the same, irrespective of whether one uses the exact endpoint exponent of 
( |2.15| ) or its variational estimate. The estimates Q2.31| - f2T33| ) for the coefficients $ differ in at 
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FIG. 1. Comparison of the two alternative mass formulae, (|2l6|) for M 2 {(3) and ( ^25|) for 
M 2 (/3), and the second order bosonization results [21] |2^] for m = 0.1. The vertical line marks 
the minimum of the curve M 2 (f3) yielding the variational estimate for f3 using 't Hooft's ansatz. 



For this value of (3 one finds M 



0.1) = 1.3969 and M 2 (m = 0.1) = 1.3913. The bosonization 



result (short-dashed horizontal line) is Mg os (m = 0.1) = 1.390±1-10~ 3 , where the error is basically 
an estimate of the unknown third order coefficient. Thus, within 't Hooft's ansatz, the alternative 
mass formula yields a somewhat better result at m = 0.1. In the next sections we will refine our 
methods so that the variational estimates for M 2 and M 2 will converge towards each other. 




least three respects from the exact values of (|2.15|) : (i) 02 and (3% depend on a, (ii) f3 2 7^ 0, 



(iii) /?3 — /3i/10. However, all these shortcomings affect at most the third order coefficient 
M 3 , and thus become negligible for small m. Only if one wants to have reliable numbers for 
M 3 , (which we will be going to produce), these effects have to be taken into account. The 
present values of $ lead to 

M 3 (a = 0) = 0.043597 , 

M 3 (a = 1) = 0.190372 , (2.34) 



which should be compared with fl2.21| ) and ( |2.24| ), respectively. The differences are large 



so that the agreement is not particularly good. These discrepancies, however, are not dis- 
turbing at this point, as we are calculating a third order effect with just one variational 
parameter. We will do better later on by enlarging the number of these parameters until we 
see satisfactory convergence of the results. 

The variational estimate for (3 can also be plugged into the alternative mass formula, 
M 2 (m,j3), ( |2.25| ), which by itself does not constitute a variational problem. We do not give 
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the analytical results for the Mi here but simply refer to Fig. 1 for a qualitative comparison 
of the two alternative formulae, and to Section |ITT| , Tables |T| and [VI|, for the actual numbers. 

As stated repeatedly, 't Hooft's ansatz ( |2.9| ) used so far has to be extended if one wants 
to accurately determine the second and third order coefficients of M 2 . This is our next issue. 



III. EXTENSION OF 'T HOOFT'S ANSATZ 



For the numerical calculations of |36| |38[ |4T| the following set of trial functions has 
been used: 

N 

<f>(x) = J2u k [x(l-x)f +k . (3.1) 

k=0 

Obviously, the term with k = and normalization w = 1 corresponds to 't Hooft's original 
ansatz ( |2.9|) . The coefficients Uk are treated as additional variational parameters, so that, 
according to the variational principle, fl3.1| ) must yield a better result than Q2.9| ). The 
question is, how big the improvement will be. To this end we will calculate the coefficients 
Mi by adding more and more basis functions to 't Hooft's original ansatz ( |2.9| ), thus enlarging 
our space of variational parameters. We will follow the two different approaches mentioned 
above, namely, (i) use the exact f3 from (|2.15 ), or, (ii) treat f3 as one of the variational 



parameters. For both approaches we will compare the calculated coefficients Mj with the 
Mi obtained from the alternative mass formula (|2.7| ) . We will continue adding basis functions 
until we see our results converge. The maximum number N of basis functions in this section 
will be five, i.e. 't Hooft's original wave function plus four corrections. 

A particular benefit of the ansatz (|3.1|) is the fact that the integrals I\ to I4 can still 



be evaluated analytically though the formulae become rather lengthy. For this reason we 
will make heavy use of the program package MAPLE in what follows. As a result we have 
analytical expressions for all the quantities we calculate. As these expressions cover pages 
and pages without being very instructive, we do not display them but rather the evaluated 
numbers. In this sense, our treatment might be called 'semi-analytical'. A major advantage, 
however, of using computer algebraic manipulations is the high calculational accuracy which 
is only limited by the maximum number of digits the machine can handle. 



A. Purely Variational Approach 

In this subsection, (3 will always be treated as a variational parameter and thus be 
obtained by minimizing the 'pion' mass squared, M 2 . In the ansatz ( |3.1|) we use up to four 
additional basis functions, so that our maximum N is four. To avoid an inflation of indices 
we rename the coefficients Uk, k — 1, . . . , 4 by a, b, c and d, respectively. Their expansion 
coefficients, defined via 

a = aim + a^m 2 , 
b = b\m + b2m 2 , 
c = Ciin + C2in 2 , 

d = dim + d 2 m 2 , (3.2) 
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are determined (recursively) by minimizing M 2 with respect to them. Let us start with the 
't Hooft model (a = 0). 



1. 't Hooft Model 

In order to save space we do not display all the values for the coefficients in ( |3.2|) . The 
best values obtained are given in the last section when we discuss the quality of the wave 
function. Here, we rather display the results for the coefficients fa of (3 (see Table |), as 
these can be compared with the exact values of Q2.15Q . 



As an important finding we note that fa remains unchanged at the exact value v^3/tt (a 
'variational invariant'), fa tends to zero, as it should. The non- vanishing of fa only affects 
the coefficient M 3 , as Mi and M 2 do not depend on fa. The convergence of fa towards 
—fa/10 seems somewhat slow; but as all Mi, i = 1,2,3, are independent of fa this has no 
observable effect. 

In Table [H] we list the expansion coefficients of the mass squared, M 2 . One notes that 
the convergence of the results for M 2 and M 3 is rather good. For M 2 we finally have a 
relative accuracy of 8 • 10~ 7 , and for M 3 of 4 • 10~ 5 . Furthermore, the coefficients are getting 
smaller if one adds more basis functions, in accordance with the variational principle. As the 
coefficient Mi is entirely determined by the 'variational invariant' fa it remains unchanged at 
27r/y / 3. This is another important result: Mi stays fixed at the standard 't Hooft-Bergknoff 
value flLJ). 

If we use the alternative mass formula (|2.7|) and evaluate it using the wave function 
calculated above we find the values for M 2 listed in Table |TTT[ 

It should be pointed out that the values of the Mj are somewhat more sensitive to the 
values of the expansion coefficients fa, because M 2 does depend on fa, and M 3 on fa and fa. 
To check the quality of our wave functions one should compare the results for Mj and Mj 
which are listed in Tables |T| and III. Once more it is gratifying to note that we are improving 



our results step by step in the variational procedure. The values of Mj converge towards 
those of Mj. The relative accuracy is approximately 10 -3 . This error is entirely due to the 
difference between the calculated and the real wave function. It is clear that the accuracy 
in the single eigenvalue M 2 (Table |I|) is higher than the one for the wave function, where 
in principle an infinite number of points have to be calculated. This represents an indirect 
proof that we are indeed improving our wave functions, and not just the eigenvalues. 



2. Schwinger Model 

For the Schwinger model (a = 1), we perform exactly the same calculations. The 
expansion coefficients of (3 are listed in Table |IV|. Although the coefficients fa, fa are 
somewhat different from those in Table p] (which they must be as they lose their dependence 
on a only in the strict limit of exact evaluation) we note the same tendency: fa converges 
to zero and fa towards — fa/10. 

The best values for the 'pion' mass squared are again provided by the variational results 
listed in Table |V|. The numerical accuracy is practically the same as for the analogous Ta- 
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TABLE I. Expansion coefficients of the end point exponent (3 for the 't Hooft model obtained 
by successively varying with respect to (3 (first line) , (3 and a (second line), etc. The first column, 
Pi, coincides with the result from the exact endpoint analysis, which in addition yields Pi = 0, 
p 3 = -ft/10. 



Ansatz Pi = \/3/tt /?2 Pz 

't Hooft 0.55132889 0.11689763 0.00490773 

a 0.55132889 0.00976951 -0.04010634 

b 0.55132889 0.00256081 -0.04889172 

c 0.55132889 0.00102086 -0.05209341 

d 0.55132889 0.00050738 -0.05343097 



TABLE II. Expansion coefficients of M 2 for the 't Hooft model obtained by successively en- 
larging the space of variational parameters. Mi is the standard 't Hooft-Bergknoff result. Note 
the good convergence towards the bottom of the table. 



Ansatz Mi = 2vr/\/3 M 2 M 3 

't Hooft 3.62759873 3.61542218 0.043597197 

a 3.62759873 3.58136872 0.061736701 

b 3.62759873 3.58107780 0.061805257 

c 3.62759873 3.58105821 0.061795547 

d 3.62759873 3.58105532 0.061793082 



TABLE III. Expansion coefficients of the alternative mass squared M 2 for the 't Hooft model 
obtained by successively enlarging the space of variational parameters. Again, the fixed value for 
Mi is the standard 't Hooft-Bergknoff result. 



Ansatz Mi = 2ir/\/3 M 2 M 3 

't Hooft 3.62759873 3.23084437 0.130791594 

a 3.62759873 3.54922830 0.066592425 

b 3.62759873 3.57265307 0.059652468 

c 3.62759873 3.57769969 0.060282959 

d 3.62759873 3.57938609 0.060854772 
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TABLE IV. Expansion coefficients of the endpoint exponent (3 for the Schwinger model ob- 
tained by successively enlarging the space of variational parameters. The first column, 0i, coincides 
with the result from the exact endpoint analysis, which in addition yields 02 = 0, = — 0i/W. 



Ansatz 0\ 02 03 

't Hooft 0.55132889 0.18250945 -0.026382222 

a 0.55132889 0.01177681 -0.040372413 

b 0.55132889 0.00317437 -0.048081756 

c 0.55132889 0.00127507 -0.051602713 

d 0.55132889 0.00063609 -0.053130398 



TABLE V. Expansion coefficients of M 2 for the Schwinger model obtained by successively 
enlarging the space of variational parameters. Again, the fixed value for Mi is the standard 
't Hooft-Bergknoff result. Note the good convergence towards the bottom of the table. 



Ansatz Mi = 2ir/V3 M 2 M 3 

't Hooft 3.62759873 3.39956798 0.19037224 

a 3.62759873 3.30906326 0.34776772 

b 3.62759873 3.30864244 0.34820193 

c 3.62759873 3.30861240 0.34820513 

d 3.62759873 3.30860791 0.34820389 



TABLE VI. Expansion coefficients of the alternative mass squared M 2 for the Schwinger model 
obtained by successively enlarging the space of variational parameters. Again, the fixed value for 
Mi is the standard 't Hooft-Bergknoff result. 



Ansatz M x = 2tv/VS M 2 M 3 

't Hooft 3.62759873 2.79913596 0.57111672 

a 3.62759873 3.27031909 0.36868190 

b 3.62759873 3.29819920 0.34869835 

c 3.62759873 3.30441759 0.34753576 

d 3.62759873 3.30651525 0.34762562 
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ble |T[ The comparison with the alternatively calculated coefficients M 2 is given in Table 
Again, everything is completely analogous to the case of the 't Hooft model (a = 0). 



B. Variational Approach using the exact (3 

In this subsection we calculate M 2 and M 2 for both values of a using the exact value 
( |2.15| ) for (3. Thus only a, b, c, d are treated as variational parameters. This is the procedure 
employed numerically in |38| . From the discussion of the preceding subsection, in particular 
the values for the (3i displayed in Table |, which differ minimally from the exact values, we 
expect that the results will be very close to those from the purely variational approach. 



1. 't Hooft Model 

This is exactly what happens as can be seen from Tables |V11| and |V1H| Comparing 
Table |VI1| with Table one finds that the coefficients Mi and M 2 of both tables coincide 
as these only depend on (3\ which is the same in both approaches. Only for M3 there are 
slight differences, due to the dependence on /3 2 . For the Mi, i = 2,3, the discrepancies are 
somewhat bigger as these coefficients do depend on /3 2 and (3% (Table |T| vs. Table [VII]] ). 
Still the consistency is quite obvious. 

Upon comparing the second columns of Tables |lj and |T|, respectively Tables [VII| and 



I Villi , one notes a funny coincidence. In the first case M 2 is slightly bigger than M 2 , and vice 
versa in the second case. If one denotes the purely variational results with a superscript V, 
and the results obtained with the exact f3 with a superscript 'e', one finds that M 2 , which 
is the same in both approaches, is given by the arithmetic mean 



Mo 



\ (M 2 e + M|) . (3.3) 



We have checked this analytically for 't Hooft's ansatz. For the higher orders this is difficult 
to do but the numerical evidence is beyond doubt. 

Comparing the third columns of Tables [TI| and |V11| , one finds that M3 < Mf. We thus 
see a slight tendency that the results of the purely variational procedure are better than 
those obtained using the exact (3. The same observation has been made by Mo and Perry 



37], in particular for large values of the fermion mass m (where the third order coefficient 
M 3 becomes important). We interpret this fact as a hint that for larger m the behaviour 
of the wave function in the intermediate region becomes more relevant compared to the 
endpoint behaviour. 



2. Schwinger Model 

For the Schwinger model, the analogous results are listed in Tables [IX| and |X|. Exactly 
the same remarks as above apply including the size of the errors and the relation between 
M 2 and Mf , M 2 . In the next section we will verify our results by purely numerical methods. 



15 



TABLE VII. Expansion coefficients of M 2 ('t Hooft model) obtained by using the exact 
endpoint exponent (3 and successively enlarging the space of variational parameters.. Again, the 
fixed value for Mi is the standard 't Hooft-Bergknoff result. 



Ansatz 


Mi = 2n/y/3 


M 2 


M 3 


't Hooft 


3.62759873 


3.61542219 


0.12513878 


a 


3.62759873 


3.58136873 


0.06230622 


b 


3.62759873 


3.58107781 


0.06184441 


c 


3.62759873 


3.58105821 


0.06180178 


d 


3.62759873 


3.58105533 


0.06179462 



TABLE VIII. Expansion coefficients of the alternatively defined mass squared M 2 ('t Hooft 
model) obtained by using the exact endpoint exponent (3 and successively enlarging the space of 
variational parameters. Again, the fixed value for Mi is the standard 't Hooft-Bergknoff result. 



Ansatz 


Mi = 2tt/V3 


M 2 


M 3 


't Hooft 


3.62759873 


4.0 


0.36275987 


a 


3.62759873 


3.61350915 


0.08659575 


b 


3.62759873 


3.58950254 


0.07106398 


c 


3.62759873 


3.58441672 


0.06603497 


d 


3.62759873 


3.58272458 


0.06406179 



TABLE IX. Expansion coefficients of M 2 (Schwinger model) obtained by using the exact 
endpoint exponent (3 and successively enlarging the space of variational parameters. Again, the 
fixed value for Mi is the standard 't Hooft-Bergknoff result. 



Ansatz 


Mi = 2t:/V3 


M 2 


M 3 


't Hooft 


3.62759873 


3.39956798 


0.38913656 


a 


3.62759873 


3.30906326 


0.34859533 


b 


3.62759873 


3.30864244 


0.34826210 


c 


3.62759873 


3.30861239 


0.34821485 


d 


3.62759873 


3.30860791 


0.34820630 
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TABLE X. Expansion coefficients of the alternatively defined mass squared M 2 (Schwinger 
model) obtained by using the exact endpoint exponent (3 and successively enlarging the space of 
variational parameters. Again, the fixed value for M\ is the standard 't Hooft-Bergknoff result. 



Ansatz M 1 = 2ir/V3 M 2 M 3 

't Hooft 3.627598730 4.0 0.3627598730 

a 3.627598730 3.347807427 0.3637551073 

b 3.627598730 3.319085690 0.3566261465 

c 3.627598730 3.312807215 0.3523062162 

d 3.627598730 3.310700593 0.3504584408 



TABLE XL Numerical results using lowest order (two-particle) LFTD approximation for 
M 2 — 1 as a function of the fermion mass m, for a = 1 (Schwinger model). 



Ansatz m = 0.0001 m = 0.0005 m = 0.001 m = 0.005 m = 0.01 m = 0.05 m = 0.1 m=0.3 m=0.5 

't Hooft 0.000362794 0.00181465 0.00363100 0.0182230 0.0366163 0.189927 0.397130 1.403677 2.705008 

a 0.000362858 0.00181465 0.00363086 0.0182207 0.0366072 0.189695 0.396181 1.394169 2.675165 

b - 0.00181461 0.00363100 0.0182210 0.0366071 0.189694 0.396176 1.394129 2.675072 

c - 0.00181459 - - 0.0366071 0.189694 0.396176 1.394125 2.675059 

d - 0.0366073 0.189694 0.396176 1.394125 2.675057 

e - - -- -- - 1.394123 2.675058 



IV. COMPARISON WITH NUMERICAL RESULTS 



For the numerical calculations the ansatz fl3.1|) has been used. The value for (3 has been 
determined numerically from ( |2.14| ). Again, up to four additional basis functions have been 
included (for larger masses even five). Let us begin with the Schwinger model (a = 1). Here 
we can use the code developed in |38|]. In Table XI we list M 2 — 1 as a function of the fermion 



mass m, calculated within two-particle LFTD approximation. The notations a,b,c,d are as 
in the preceding section, the letter 'e' denotes the inclusion of a fifth basis function. 

One main difference in comparison with the computer algebraic treatment of Section [Tl I 
is the numerical inaccuracy for small m. This is a general disease of numerical treatments, 
and also shared e.g. by the lattice [|T3J or DLCQ approach P3fl . In these approaches, 
however, the numerical errors are typically much larger than ours (see the next section for 
an explicit comparison). Due to the small-m instability, our numerical calculation does not 
converge within an arbitrary number of digits. As soon as a calculated value for M 2 — 1 
becomes bigger as the preceding one (a numerical violation of the variational principle), 
we terminate the procedure and pick the smaller value as our final result. The difference 
between these last two values can be used as an estimate of the numerical inaccuracy. The 
errors will be further discussed below and in the next section when we compare our results 
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with related work (see Tables |XV| , |XV1| ). 

It should also be pointed out that the numerical analysis is conceptually very different 
from the computer algebraic treatment. As one cannot perform a Taylor expansion nu- 
merically, the expansion coefficients of M 2 have to be obtained by fitting polynomials to 
M 2 (m). Clearly, this is an additional source of errors, and one expects the results to be less 
accurate than those of the preceding sections. A cubic fit to the optimum values for M 2 — 1 
in Table |XI] yields the following expansion coefficients of M 2 , 

M 1 = 3.62609 , (4.1a) 
M 2 = 3.33607 , (4.1b) 

M 3 = 0.22396 , (4.1c) 

which should be compared with the last line of Table [V|. To estimate the accuracy of these 
values we also show the results from a quartic fit, 

Mi = 3.62755 , (4.2a) 
M 2 = 3.31029 , (4.2b) 

M 3 = 0.32984 . (4.2c) 

Comparing ( |4.1|) and (|4.2j ) one finds that the stability of the fits is not too impressive in 
view of the accuracy we would like to achieve. In particular the third order coefficient is 
numerically difficult to determine. We estimate the relative accuracy as being 10~ 3 for Mi, 
10~ 2 for M 2 and 10 _1 for M3. This is also confirmed by comparing with Table [V]. 

For the Schwinger model it is possible to check the influence of higher particle sectors on 



the 'pion' mass squared by using the machinery developed in |38[ for the wave functions of 
the higher Fock components. These are the amplitudes of finding not only two, but four, six, 
... (anti-)fermions in the 'pion'. In Table |XII| , we list the best values for M 2 — 1 including up 
to six-body wave functions. We also show the 2-, 4-, and 6-body content of the total wave 



function. It is known that the 'pion' is entirely 2-particle in the chiral limit [j37 |. For small 
mass, one therefore expects only small contributions from the higher Fock sectors. This is 
confirmed by the numerical results. Astonishingly, this feature persists up to values of at 
least m = 0.5 for the fermion mass. This fact is in agreement with the observation of Mo 
and Perry that the four-particle component of the wave function is less than 0.4 % for all 
values of the fermion mass [^7f| . A similar result has been found in the DLCQ calculations of 
33| . In addition we note that there seems to be some interesting kind of hierarchy between 



the relative strengths of the contributions from different particle sectors. If we denote the 
2/c-particle amplitude in the wave function by f 2 k, we find that \f 2 \ 2 ^> I/4I 2 \f6\ 2 , the 
individual proportions being several orders of magnitude (see Table |XHj ) . Here, we explicitly 



see the magic of light-front field theory at work: high Fock components in bound states tend 
to be largely suppressed, at variance with the situation encountered in field theory quantized 
the standard way. 
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TABLE XII. Numerical Results for M 2 — 1 function of the fermion mass m, for a = 1 
(Schwinger model) including 2-, 4- and 6-body wave functions. The contributions of the different 
Fock sectors to the total wave function squared are given in percent. 



m = 0.0001 m = 0.0005 m = 0.001 m = 0.005 m= 0.01 m = 0.05 m = 0.1 m = 0.5 

M 2 - 1 0.000362793 0.00181481 0.00363055 0.0182179 0.0365968 0.189468 0.395400 2.66787 

% 2-body 100.00 100.00 100.00 99.999990 99.999962 99.999213 99.997514 99.986816 

% 4-body - 0.000010 0.000038 0.000782 0.002474 0.013171 

% 6-body - - 0.000005 0.000013 0.000013 



From the point of view of the variational principle, the results shown in Table [XI 1| are a 
bit better {i.e. smaller) than those of Table |X| (apart from the value for m = 0.0005), but 
the improvement is rather small. A cubic fit yields 



M = 3.62667 



(4.3a) 



Mo = 3.23696 



(4.3b) 



and a quartic fit 



M, = 0.36235 



Mi = 3.62747 



(4.3c) 



(4.4a) 



Mo = 3.20864 



(4.4b) 



M, = 0.60365 



(4.4c) 



Compared with (4.1) and (|4.2|) we do not find any absolute improvement in Mi that 
could be distinguished from the numerical inaccuracy. For the coefficient M2, which in 
bosonization schemes is 3.3874, we even get the wrong tendency: it becomes smaller upon 
including higher Fock states. The inaccuracy for M 3 is so large that this coefficient is only 
determined in its order of magnitude. 

Altogether, we arrive at the very important conclusion that the inclusion of higher 
particle sectors in the light-front bound state equation of the Schwinger model does not 
diminish the discrepancy between the results obtained via the 't Hooft-Bergknoff method 
and those from bosonization techniques. 

Let us move to the 't Hooft model. If we put a = in the Schwinger model code above 
(Code I) we find the 'pion' mass squared for the 't Hooft model. The best values within 
the 2-particle sector are listed in Table Xill . Note that higher particle sectors are strictly 
suppressed in the large- N c limit. Thus, there is no point in calculating these, unless one is 
interested in 1/N C corrections to the 't Hooft model results |2|| |34| [ |35[ . 



In view of the restricted number of data points only a quadratic fit makes sense which 
yields 
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TABLE XIII. M 2 as a function of the fermion mass m, for a = ('t Hooft model), obtained 
with Code I 



m 


0.0001 


0.001 


0.01 


0.1 


M 2 


0.000362795 


0.00363109 


0.0366342 


0.398634 



TABLE XIV. M 2 as a function of the fermion mass m, for a = ('t Hooft model), obtained 
with Code II 



771 


0.0001 


0.001 


0.01 


0.1 


M 2 


0.000362795 


0.00363118 


0.0366341 


0.398634 



M l = 3.62754 , (4.5a) 

M 2 = 3.58796 , (4.5b) 

which is consistent with the results of Section [Til] (see Table |IT]). A second code (Code II), 
which was independently developed for the 't Hooft model |35| yields the results displayed 
in Table [X1V| . From a quadratic fit we obtain 

M 1 = 3.62754 , (4.6a) 

M 2 = 3.58806 . (4.6b) 

Comparing with (f4.5|) we see that both codes yield the same results within the numerical 
accuracy. This is reassuring, since, as already stated, the codes were developed independently 
from each other. 

A cubic fit to a total of 30 data points produced with Code II yields 

M l = 3.62758 , (4.7a) 



M 2 = 3.58260 , (4.7b) 

M 3 = 0.06450 , (4.7c) 

in fair agreement with the variationally obtained values of Table O. To check the stability 
of this fit we can compare with the extension to fourth order, 

M 1 = 3.62756 , (4.8a) 



M 2 = 3.58314 
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(4.8b) 



M 3 = 0.06280 , (4.8c) 

M 4 = 0.001232 . (4.8d) 

This leads to an estimate of roughly 10~ 6 , 10 -4 and 10~ 2 for the (absolute) numerical error 
in the first, second and third order coefficient, respectively. Furthermore it is gratifying to 
note that the fourth order coefficient, M4, is numerically small. 

As is obvious from the discussion above, the numerical errors for the 't Hooft model are 
much smaller than those for the Schwinger model. The basic reason for this is the Schwinger 
model anomaly. Note that in the bound state equation the anomaly factor a multiplies an 
integral over the wave function. To evaluate this, one needs the wave function as a whole. In 
the 't Hooft model, on the other hand, this term is absent and the eigenvalues are dominated 
entirely by the endpoint behaviour of the wave function which is known exactly. This makes 



the numerical errors much smaller (see also Tables [XV], |XVI|) . 



V. DISCUSSION AND CONCLUSION 

In the preceding sections we have calculated the expansion coefficients Mj in the series 

M 2 = a + M x m + M 2 m 2 + M 3 m 3 (5.1) 

for the 'pion' mass squared of both the 't Hooft (a = 0) and the Schwinger model (a = 1). 
In order for the expansion ( |5.1| ) to make sense we have considered only small masses m < 1 
i.e. m <C fi in the original units. For the 't Hooft model we have /i 2 , = g 2 Nc/2ir with 
Nc — > 00, g 2 Nc fixed. Thus, we are working in the weak coupling phase where the limit 
No — > 00 (g — ► 0), is taken before the limit m — > 0, or, equivalently, such that one always 
has m ^> g ~ l/y/N c -^00 jEJ. For the Schwinger model, /1q = e 2 /Tr, so that small 
fermion mass corresponds to strong coupling. 

We have used analytical, computer algebraic and numerical methods of different accu- 
racy. Nonetheless, the overall picture is intrinsically consistent. Our main findings are the 
following: 

(i) The first order coefficient, Mi, in the expansion of the 'pion' mass squared is independent 
of the anomaly a, i.e. the same in both the 't Hooft and Schwinger model. This confirms 
the results of 't Hooft § § and Bergknoff |0. 

(ii) The 't Hooft-Bergknoff value, M\ = 27r/v3, is a 'variational invariant': it does not get 
altered by extending the space of variational parameters. This has been checked analytically 
and numerically. 

(iii) In the Schwinger model, the 't Hooft-Bergknoff value does not change upon inclusion 
of higher particle sectors. Only the second and third order coefficients, M 2 and M 3 , are 
affected. 

(iv) Thus, for the Schwinger model, there remains a few percent discrepancy in the coeffi- 
cients Mi and M 2 compared to bosonization results. 

(v) The variational calculation yields the 'pion' wave function to a high accuracy. The 
endpoint behaviour is reproduced exactly (in leading order in m). The behaviour in the 
intermediate region, < x < 1, gets improved as can be seen from Fig. 3 below. 
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TABLE XV. The expansion coefficients Mj for the 't Hooft model (a = 0). The errors of 
our results obtained within 2-particle Tamm-Dancoff (2PTD) approximation are estimated by 
comparing the last two lines in Table |J (for the variational method) and the different polynomial 
fits ( [4.5| - [4.8|) (for the numerical results). The numerical inaccuracies given in the last line were 
estimated by us using the polynomial fit method. 



Mi 



Mo 



M, 



variational 2PTD 
numerical 2PTD 
't Hooft [§ 
Burkardt |26l 
Li |2l 



Li et al. @ 



2n/^/S = 3.627599 
3.62758 ±2 • 1(T 5 
27T/V3 = 3.627599 
2ir/V3 = 3.627599 
3.64 ± 0.06 
3.64 ±0.03 



3.581055 ±3 • 10~ 6 
3.5829 ±3- 10~ 4 

3.5812 

3.60 ±0.06 



0.061793 ±3 • 10~ 6 
0.064 ±1 • 10~ 3 



0.04 ±0.04 



In Table [XV| we summarize our optimum final results for the 't Hooft model (a = 0) 
and compare with results that have been obtained previously. 

Upon inspection of Table |XV| one finds a very good overall consistency of the results. 
All given values are in good agreement within error bars. The scale is set by the variational 
results which are the most accurate ones. They are matched by our numerical results as 
well as by previous analytical and numerical calculations. 

As far as the latter are concerned, a few remarks are in order. The coefficients Mj 
displayed in the last line of Table [XV] are obtained by performing polynomial fits to the 
values calculated numerically in ^3]] (which basically agree with those from 't Hooft 's original 
calculation displayed in Fig. 5 of [0]). The averaged result of the fits is shown including an 
estimate of the errors. It should also be stressed that the data in [13] are mainly obtained for 
m > 0.5 whereas the bulk of our values for the Mj are obtained for m < 0.1. Nevertheless, 
the agreement with our results is satisfactory within error bars. 

From Table [XV| it is obvious that all numerical results favour the first order 't Hooft- 
Bergknoff value of 2n/ \^3. This agreement is particularly gratifying for the data of [43 which 
were not obtained via LC techniques but within ordinary quantization and even within a 
different gauge, namely axial gauge. The agreement is therefore highly non-trivial [fii]] . 

For the Schwinger model, (a = 1), we summarize our findings in Table [XV1| . For 
comparison we have also listed analytical results obtained via bosonization techniques and 
data which we extracted from polynomial fits to numerical LFTD []3"7|] , DLCQ []33J and lattice 
results [45). The contents of Table |XVI| are graphically displayed in Fig. 2, where we have 



chosen rescaled units as in 22 



As already stated, due to the anomaly the numerical errors for the Schwinger model 
are at least an order of magnitude larger than for the 't Hooft model. In addition, the 
2% discrepancy between our LC results and the analytic bosonization results does not get 
resolved. However, we have shown that the discrepancy is not due to (i) inaccuracies in the 
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FIG. 2. The rescaled 'pion' mass M = M/y/l + m 2 as a function of the fermion mass 
m. The short-dashed curve represents a 'phenomenological' parametrization for M with 
M 2 = l+Mim+4m 2 , which becomes exact for very small and very large m and thus smoothly inter- 
polates between the strong and weak coupling regions. The long-dashed curve is our second-order 
and the solid curve our third-order result. As expected, mass perturbation theory breaks down as 
m becomes of order 1. The crosses are the lattice results 



J , the diamonds the lattice results |46[ , 
which go down to comparatively small masses, however with large errors. The circles are LFTD 



results |p7 |, whereas the squares |33| and triangles |47|, as quoted in p2 |, are DLCQ results. The 
values corresponding to the triangles are not included in Table XVI| . The 2% discrepancy in Mi is 
invisible within the resolution of the figure. 
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0.8 
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m 100. 
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TABLE XVI. The expansion coefficients Mj for the Schwinger model (a = 1). The errors for 
the variational 2-particle Tamm-Dancoff (2PTD) approximation are estimated by comparing the 
last two lines in Table [V|. The errors of our numerical calculations are obtained from comparing 



polynomial fits of different order to the numerical results of Tables |XI| and [XI]] . Errors of other 
results are given where they could be estimated analogously. 



Mi 



M 2 



M 3 



variational 2PTD 
numerical 2PTD 
numerical 6PTD 
numerical 4PTD 



DLCQ || 



lattice m 
Adam |2l| 
Fields et al. fl22l 



2it/^3 = 3.627599 
3.6268 ±8 • 10" 4 
3.6267 ±4 • 10" 4 
3.62 ±0.07 
3.7 ±0.2 
3.5 ±0.2 
2e 7 = 3.562146 
2e 7 = 3.562146 



3.308608 ±4-10- 
3.32 ±0.02 
3.22 ±0.02 
3.2 ±0.3 
3.5 ±0.3 
3.7 

3.3874 
3.387399 



0.348204 ±1 • 10~ 6 
0.28 ±0.06 
0.5 ±0.1 
0.3 ±0.2 

0.02 



wave function or (ii) neglect of contributions from the next higher Fock sectors (four and 
six particles). It should be mentioned that other numerical methods (DLCQ, lattice) are by 
far too inaccurate to distinguish between the 't Hooft-Bergknoff and the bosonization value. 
We will come back to these issues in a moment. 

As is well known, p2| |37] [0, the DLCQ data are comparatively inaccurate for small 
fermion mass m (i.e. large coupling). This is due to the fact that the dominating feature of 
the LC wave function, its endpoint behaviour, is not very accurately reproduced using an 
equally-spaced momentum grid. The poor convergence of DLCQ for small m has recently 
been overcome by incorporating the exact endpoint behaviour [plHf . 

The lattice data generally suffer form the same disease of having rather large error bars 
for small m. The chiral limit, m = 0, can only be reached via extrapolation in a very 
inaccurate manner unless one uses the Schwinger result M 2 (m = 0) = 1 as a bias [45]. A 



representative collection of lattice results for Mi is given in Table |X.V11| . For the convenience 
of the reader we also list the condensate (0|^V|0) = ~~ Mi/ An. 

The values for the Schwinger model condensate in Table |XVII| should be compared 
with our results. If we assume that formulae (O) and (L7) are directly applicable to our 
calculation we get 

(O|<0V|O) = - — = ^7= = -0.28868 , (5.2) 

4vr 2y/3 

1 / Mo 

(<wio> = -^yxT5 = -°- 28015 ' ( 5 - 3 ) 

while the standard result is 
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TABLE XVII. Lattice results for the first order coefficient Mi and the (negative of the) 
condensate, —{0\^ip\0} = Mi /Air. The results are quoted in chronological order. They should be 
compared with the bosonization results, Mi = 3.562, —{0\^ip\0) = 0.283 and the 't Hooft-Bergknoff 
values, Mi = 3.628, -(O|^0|O) = 0.289. 





BKS76 a 


CKS76 b 


CH80 C 


MPR8I d 


HKC82 C 


CK86 f 


Mi 


3.42 


3.644 


3.48 


2.97 


3.33 


3.77 


-<o|V^|o> 


0.27 


0.290 


0.28 


0.24 


0.26 


0.30 



a Banks, Kogut, Susskind [12] 

b Carroll, Kogut, Sinclair, Susskind |4 

c Crewther, Hamer [45 

d Marinari, Parisi, Rebbi [15] 

c Hamer, Kogut, Crewther, Mazzolini 

f Carson, Kenway [46] 



(0|^|0) = -— = -0.28347 . (5.4) 

As already mentioned in the introduction, the few percent discrepancies in Mi and M 2 
immediately affect the condensate. Note that the condensate value obtained from the second 
order coefficient is closer to the bosonization result, the relative error being 1.2 % compared 
to the 1.8 % at first order. Furthermore, the exact result ( |5.4|) lies between our first and 
second order values, ( |5.2[ ) and ( |5.3| ), so that their mean value, (0|^ , ?/'|0) = —0.28442, is 
much closer to ( |5.4j ), the error being only 0.3%. A similar reduction of errors is actually 
at work if one considers the 'pion' mass-squared of the Schwinger model as a function of 
m. The first order overshoots the bosonization value while the second order contribution 
is too small. Adding both one gets closer to the exact result, at least if the mass m is not 
too tiny. For the value of Fig. 1, m = 0.1, one finds to second order M 2 (0.1) = 1.39585, 
while Mg os (0.1) = 1.39009. The relative difference is only 0.4%. Possible resolutions of the 
discrepancies between (|5.2|), ( |5.3| ) and fl5T4] ) will be discussed at the end of this section. 

Apart from the 'pion' mass squared, the eigenvalue in the 't Hooft-Bergknoff equation, 
we have also calculated the associated eigenfunction, the light-cone wave function, with high 
accuracy. The latter has been obtained with high accuracy as can be seen from the good 
convergence of the alternatively defined mass squared coefficients, M; towards the variational 
estimates, M; (see Section Q). If we denote (f>o(x) = x{l — x), our most accurate variational 
ansatz for the wave function can be written as 

[fa a, b, c, d] = $ + a 4 +1 + b <^ +2 + c <^ +3 + d <^ +4 . (5.5) 

According to our mass perturbation theory, each of the variational parameters is expanded 
in powers of the fermion mass m, the coefficients being denoted fa, fa, fa, a i, a 2, etc. From 
't Hooft's endpoint analysis, (3 is known exactly. It turns out to be independent of the 
anomaly (a) and thus is the same for the 't Hooft and Schwinger model. In Table |XVI i I 



we compare the exact expansion coefficients of j3 with their variational estimates. The best 
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TABLE XVIII. Comparison of the exact expansion coefficients of the endpoint exponent 
with their best variational estimates for both the 't Hooft (a = 0) and the Schwinger model (a = 
!)• 





01 


02 


03 


exact 


^/3/vr = 0.55132890 





-P 1 /10 = -0.05513289 


variational (a = 


0) V3/ir = 0.55132890 


0.00051 


-0.053 


variational (a = 


1) v^/vr = 0.55132890 


0.00063 


-0.053 


TABLE XIX. 
the 't Hooft (a = 


Best estimates for the variational parameters 
= 0) and Schwinger model (a = 1). 


in the LC wave function for both 




a 2 b\ b 2 


c\ 


c 2 di d 2 


a = 0.92 


1.19 -2.2 -3.0 


7.0 


9.2 -9.9 -11.8 


a = 1 1.42 


1.03 -2.7 -3.1 


8.8 


9.6 -12.5 -12.1 



estimates for the other variational parameters are listed in Table |XIX| for both the 't Hooft 
and Schwinger model. 

The errors in the coefficients are comparatively large and growing from the left to the 
right in Table [X1X|. Good numerical convergence is evident for the coefficients of a. For 



the other coefficients our method does not provide enough iteration steps to make conver- 
gence explicit. However, the sensitivity of the mass expansion coefficients Mj and the wave 
functions on the parameters b, c, d is very weak (see Fig. 3). 

Let us finally investigate possible sources for the few-percent discrepancies in the 
Schwinger model results like for the condensate values, ( |5.2| ) and ( |5.3| ) vs. (|5.4|). First of all, 
it has to be reemphasized that the light-cone calculations are conceptually rather different 
from the usual treatment based on bosonization within standard equal-time quantization. 
We are studying the lowest lying meson as a relativistic fermion-anti-fermion bound state. 
The fermionic degrees of freedom are explicitly taken into account as they define our Fock 
basis. In the bosonized theory, which is a sine-Gordon model with the Lagrangian || || 

I i 

C = -<9 M 0<9 M — -/Uq0 2 — cm/i cos V47T0 , (5.6) 

the 'pion' is the elementary particle described by the scalar field <fi. One calculates its 
mass by perturbation theory in the 'coupling' cm/io, c = e 1 /2ir. In the bosonized theory, 
this 'coupling' is actually dependent on the normal-ordering scale used to renormalize the 
theory. The standard choice is the most natural one, namely the Schwinger boson mass /io, 
which explicitly appears in the Lagrangian ( |5.6|) . It is, however, not very hard to determine 
the general scale dependence of quantities calculated within perturbation theory. All one 
needs is Coleman's re-normal-ordering prescription |5(J which for the case at hand is 
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FIG. 3. The light-cone wave function of the 't Hooft model 'pion' for m = 0.1. The solid 
curve represents the result from 't Hooft 's original ansatz, the dashed curve our best result (with 
maximum number of variational parameters). At the given resolution, however, the curves of all 
extensions of 't Hooft's ansatz (a, b, c, d) lie on top of each other. 




iV/io cos 



AlT(f) 



V AT 

— 7V„ cos 

1*0 



4-7T0 



(5.7) 



Here, iV M0 and iV M denote normal-ordering with respect to /xq and fi, respectively. From 
one derives the following identity for the boson mass squared (setting /j = 1), 



M 



1 + Mi m + M 2 m 2 + M 3 m 3 + . . . 

1 + M[ nm + M' 2 /i 2 m 2 + M' z fx 3 m 3 + ... . 



(5.8) 



In this expression, the Mj denote the coefficients calculated with the standard normal- 
ordering scale /i , where all tadpoles can be set to zero, while the M[ are calculated with 
normal-ordering scale /i, where one has non-vanishing tadpole contributions. In order that 
the 'pion' mass be scale independent one has to have that (upon reintroducing /i ) 



M' 

n 



(5.9) 



We have checked this identity explicitly in the first two non-trivial orders of mass per- 
turbation theory with the Lagrangian (|5.6|). It is also interesting to note that the scale 
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dependence of the coefficients might equivalently be associated with a multiplicative rescal- 
ing of the fermion mass m by defining m' = /xm and using ml as the expansion parameter 

Now, the upshot of all this is the following. As we are working entirely in the fermionic 
representation, we do not a priori know, to which normal-ordering scale of the bosonized 
theory our results correspond. If we assume that it is not the Schwinger mass /Xo but a 
different scale /x, we can match the first order coefficients by choosing 

" = m "° = ' (5 - 10) 

where M[ now denotes our light-front result. Quite a similar point of view has been taken by 
Burkardt in his paper on light-cone quantization of the sine-Gordon model j5l|. However, 



as is obvious from ([5.9|), the change of scale ( |5.10| ) propagates systematically through all 
coefficients. In particular, if M[ is larger than the corresponding bosonization result, the 
same has to be true for all M[. Our coefficient M' 2 , however, already fails this requirement. 
We thus conclude that a possible scale dependence alone cannot explain the discrepancy. 

An additional source of error could still be contributions from higher Fock sectors. We 
have shown numerically that the first few low orders have very little impact on the expansion 
coefficients Mi. The contribution from each Fock sector might thus be exponentially small. 
The whole series, however, if it could be summed up, might add something seizable. Let us 
assume, for example, that the contribution of higher Fock sectors to M2 has the following 
form[| 

00 

5M 2 = 5> 2n (m), (5.11) 

n=2 

where 2n is the label of the 2n-particle sector. If we had C2 n (m) = e _2nm , summation of 
the series would yield 5M 2 = I /2m + 0(1) and thus effectively a contribution to the first 
order in m. As we are not able to calculate the real high order behaviour of the series in 
( |5.11| ), we cannot, at the moment, make any definite conclusion about a possible relevance 
of such summation effects. All we can say is, that if such effects are present, they must be 
numerically small, i.e. at the few percent level. 

As a first step towards getting some better control of higher particle sectors one would 
like to develop some tools to estimate their contributions analytically. This amounts to 
determine the coefficients C2 n (m), at least for small n. So far, higher Fock states have only 
been included numerically. 

It might also be interesting to calculate excited states. For the Schwinger model, this 
again requires the inclusion of higher Fock components as pointed out by Mo and Perry 
P7J . For the 't Hooft model, on the other hand, one can stay in the two-particle sector. 
One might then try to modify the variational procedure, basically by taking into account 
the growing number of nodes, with the aim to match the low energy spectrum with the high 
energy behaviour derived by 't Hooft. Work in these directions is underway. 



2 We thank M. Burkardt for suggesting the following example. 
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